Quantum simulation of the single-particle Schrodinger equation 
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The working of a quantum computer is described in the concrete example of a quantum simulator 
of the single-particle Schrodinger equation. We show that a register of 6 — 10 qubits is sufficient 
to realize a useful quantum simulator capable of solving in an efficient way standard quantum 
mechanical problems. 



I. INTRODUCTION 



During the last years quantum information 1 ' 2 has be- 
come one of the new hot topic field in physics, with the 
potential to revolutionize many areas of science and tech- 
nology. Quantum information replaces the laws of classi- 
cal physics, applied to computation and communication, 
with the more fundamental laws of quantum mechanics. 
This becomes increasingly important due to technological 
progress reaching smaller and smaller scales. Extrapolat- 
ing the miniaturization process of integrated circuits one 
would estimate that we shall reach the atomic size for 
storing a single bit of information around the year 2020 
(at the present time the typical size of circuit compo- 
nents is of the order of 100 nanometers). At that point, 
quantum effects will become unavoidably dominant. 

The final aim of quantum computation is to build a 
machine based on quantum logic, that is, a machine that 
can process the information and perform logic operations 
in agreement with the laws of quantum mechanics. In 
addition to its fundamental interest, a large scale quan- 
tum computer, if constructed, would advance computing 
power beyond the capabilities of classical computation. 

The purpose of the present paper is to illustrate the 
main features of a quantum computer in the concrete 
example of a quantum simulator of the single-particle 
Schrodinger equation. We wish to show that quantum 
logic, already with a small number of qubits, allows the 
efficient simulation of basic quantum mechanical prob- 
lems. A quantum simulator with 6—10 qubits could 
produce a real picture book of quantum mechanics. We 
believe that snapshots from this book will help develop a 
physical intuition of the power of quantum computation. 

Our paper adds to other existing pedagogical presenta- 
tions of various apects of quantum computation, includ- 
ing Shor's quantum algorithm for integer factorization 3 , 
Grover's quantum search algorithm 4,5 , quantum infor- 
mation processing by means of nuclear magnetic reso- 
nance techniques 6 , using linear optics 7 , with cold atoms 
and ions 8 , or with superconducting circuits 9 , quan- 
tum measurements 10 , decoherence 11 , quantum error- qubit, and X)fe=o 



II. QUANTUM LOGIC 

The elementary unit of quantum information is the 
qubit (the quantum counterpart of the classical bit) and 
a quantum computer may be viewed as a many-qubit 
system. Physically, a qubit is a two-level system, like 
the two spin states of a spin-i particle, the polarization 
states of a single photon or two states of an atom. 

A classical bit is a system that can exist in two dis- 
tinct states, which are used to represent and 1, that 
is, a single binary digit. The only possible operations 
(classical gates) in such a system are the identity (0 — ► 0, 
1 — ► 1) and NOT (0 — ► 1, 1 — > 0). In contrast, a quantum 
bit (qubit) is a two-level quantum system, described by 
a two-dimensional complex Hilbert space. In this space, 
one may choose a pair of normalized and mutually or- 
thogonal quantum states, called |0) and |1), to represent 
the values and 1 of a classical bit. These two states form 
a computational basis. From the superposition principle, 
any state of the qubit may be written as 



a\Q)+(5\l), 



(1) 



where the amplitudes a and (3 are complex numbers, con- 
strained by the normalization condition \a\ 2 + \(3\ 2 = 1. 
Moreover, state vectors are defined up to a global phase 
of no physical significance. Therefore, \yb) only de- 
pends on two real parameters and we may write \ip) — 



cos 



|0> 



J<t> 



sin with < 6 < tt and < cj> < 2tt. 



A collection of n qubits is known as a quantum register 
of size n. Its wave function resides in a 2 n -dimensional 
complex Hilbert space. While the state of an n-bit classi- 
cal computer is described in binary notation by an integer 

fee {0,1,..., 2™- 1}, 



k — k n —i 2 



h 2 + k , 



(2) 



with fco, fci, . . . , fc n -i S {0, 1} binary digits, the state of 
an n-qubit quantum computer is 



(3) 



k=0 



where \k) = |fe„-i) 



correction and fault-tolerant quantum computation 
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|Cfc| 

hand notation \k n -i) 



\ki)\ko), with \kj) state of the j-th 
2 = 1. Note that we use the short- 
■ •• |fci)|feo) for the tensor product 
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\k n -i) ® • • • ® |fei) ® |fco). Taking into account the nor- 
malization condition and the fact that a quantum state 
is only defined up to an overall phase, the state of an 
n-qubit quantum computer is determined by 2 x 2™ — 2 
independent real parameters. 

The superposition principle is clearly visible in Eq. (3): 
while n classical bits can store only a single integer k, the 
n-qubit quantum register can be prepared in the corre- 
sponding state of the computational basis, but also 
in a superposition. We stress that the number of states 
of the computational basis in this superposition can be 
as large as 2 n , which grows exponentially with the num- 
ber of qubits. The superposition principle opens up new 
possibilities for efficient computation. When we perform 
a computation on a classical computer, different inputs 
require separate runs. In contrast, a quantum computer 
can perform a computation for exponentially many in- 
puts on a single run. This huge parallelism is the basis 
of the power of quantum computation. 

We stress that the superposition principle is not a 
uniquely quantum feature. Indeed, classical waves satis- 
fying the superposition principle do exist. For instance, 
we may consider the wave equation for a vibrating string 
with fixed endpoints. It is therefore also important to 
point out the role of entanglement for the power of quan- 
tum computation, as compared to any classical compu- 
tation. Entanglement is arguably the most spectacular 
and counter-intuitive manifestation of quantum mechan- 
ics, observed in composite quantum systems: it signifies 
the existence of non-local correlations between measure- 
ments performed on particles who interacted in the past 
but now are located arbitrarily far away. Mathemati- 
cally, we say that a two-particle state \tp) is entangled, or 
non-separable, if it cannot be written as a simple product 
|fci)|fc2) of two states describing the first and the second 
subsystem, respectively, but only as a superposition of 
such states: — J2 kl k2 Ch 1 k 2 \ki)\k2)- For instance, the 
(Bell) state -^|(|00) + |11)) is entangled, while the state 
-^|(|00) + 1 1 0} ) is separable (i.e., not entangled), since 
we can write it in the form -^|(|0) + |1))|0). 

There is no entanglement in classical physics. There- 
fore, in order to represent the superposition of 2™ levels 
by means of classical waves, these levels must belong to 
the same system. Indeed, classical states of separate sys- 
tems can never be superposed. The overall state space 
of a composite classical system is the Cartesian product 
of the individual state spaces of the subsystems, while 
in quantum mechanics it is the tensor product. Thus, 
to represent the generic n-qubit state of (3) by classical 
waves we need a single system with 2™ levels. If A is the 
typical energy separation between two consecutive lev- 
els, the amount of energy required for this computation 
is given by A2™. Hence, the amount of physical resources 
needed for the computation grows exponentially with n. 
In contrast, due to entanglement, in quantum physics 
a general superposition of 2™ levels may be represented 
by means of n qubits. Thus, the amount of physical re- 



sources (energy) grows only linearly with n. 

To implement a quantum computation, we must be 
able to control the evolution in time of the many-qubit 
state describing the quantum computer. As far as the 
coupling to the environment may be neglected, this evo- 
lution is unitary and governed by the Schrodinger equa- 
tion. It is well known that a small set of elementary 
logic gates allows the implementation of any complex 
computation on a classical computer. This is very im- 
portant: it means that, when we change the problem, 
we do not need to modify our computer hardware. For- 
tunately, the same property remains valid for a quan- 
tum computer. It turns out that, in the quantum circuit 
model, each unitary transformation acting on a many- 
qubit system can be decomposed into (unitary) quantum 
gates acting on a single qubit and a suitable (unitary) 
quantum gate acting on two qubits. Any unitary op- 
eration on a single qubit can be constructed using only 
Hadamard and phase-shift gates. The Hadamard gate is 
defined as follows: it turns |0) into (|0) + |l))/\/2 and 
|1) into (|0) - |1))/V2. The phase-shift gate (of phase 
S) turns |0) into |0) and |1) into e t(5 |l). Since global 
phases have no physical meaning, the states of the com- 
putational basis, |0) and |1), arc unchanged. However, 
a generic single-qubit state a\0) + /3\1) is mapped into 
c*|0) + (ie l5 \l). Since relative phases are observable, the 
state of the qubit has been changed by the application of 
the phase-shift gate. We can decompose a generic uni- 
tary transformation acting on a many-qubit state into 
a sequence of Hadamard, phase-shift and controlled-not 
(CNOT) gates, where CNOT is a two-qubit gate, defined 
as follows: it turns |00) into |00), |01) into |01), 1 10) into 
|11) and 1 11) into 1 10) . As in the classical XOR gate, the 
CNOT gate flips the state of the second (target) qubit if 
the first (control) qubit is in the state |1) and does noth- 
ing if the first qubit is in the state |0). Of course, the 
CNOT gate, in contrast to the classical XOR gate, can 
also be applied to any superposition of the computational 
basis states. 

The decomposition of a generic unitary transformation 
of an n-qubit system into elementary quantum gates is in 
general inefficient, that is, it requires a number of gates 
exponentially large in n (more precisely, 0(n 2 4") quan- 
tum gates). However, there are unitary transformations 
that can be computed efficiently in the quantum circuit 
model, namely by means of a number of elementary gates 
polynomial in n. This is the case in many computational 
problems of interest. A very important example is given 
by the quantum Fourier transform, mapping a generic 
n-qubit state ^2k=o f(k)\k) into J^f^ 1 f{l)\l), where 
N = 2™ and the vector {/(0), ...,f(N — 1)} is the dis- 
crete Fourier transform of the vector {/(0), f(N — 1)}, 
that is, f(l) = ^Efjo 1 e 2mkl / N f(k). It can be 
shown that this transformation can be efficiently im- 
plemented in 0(n 2 = (log 2 N) 2 ) elementary quantum 
gates, whereas the best known classical algorithm to sim- 
ulate the Fourier transform, the fast Fourier transform, 
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FIG. 1: A quantum circuit implementing the quantum Fourier 
transform (except for the fact that the order of qubits in the 
output is reversed, so that one has simply to relabel qubits 
in the output state). As usual in the graphical representation 
of quantum circuits, each line corresponds to a qubit and any 
sequence of logic gates must be read from the left (input) 
to the right (output). From bottom to top, qubits run from 
the least significant (ko, according to binary notation (2)) to 
the most significant (fc n _i). Here a qubit is said to be more 
significant than another if its flip gives a larger variation in the 
integer number coded by the state of the n qubits. A square 
with H written inside stands for the Hadamard gate, while 
for controlled- phase shift gates CPHASE(|f-) gates we draw 
a circle on the control qubit and a square with Rk written 
inside on the target qubit. Note that to implement the inverse 
Fourier transform it is sufficient to run the circuit in this figure 
from right to left and with R j. instead of R k . This is due to the 
facts that the Hadamard gate is self-inverse and R{ = R^ 1 . 

requires 0(N\og 2 N) elementary operations. A quan- 
tum circuit computing the quantum Fourier transform is 
shown in Fig. 1: it requires n Hadamard and n(n — l)/2 
controlled phase-shift gates. By definition, the controlled 
phase-shift gate CPHASE(<5) applies a phase shift 5 to the 
target qubit only when the control qubit is in the state 
|1): it turns |00) into |00), |01) into |01), 1 10) into |10) 
and 1 11) into e J(5 |ll). The quantum Fourier transform 
is an essential subroutine in many quantum algorithms, 
including the quantum simulator described in this paper. 

III. QUANTUM ALGORITHMS 

As we have seen, the power of quantum computation is 
due to the inherent quantum parallelism associated with 
the superposition principle. In simple terms, a quantum 
computer can process a large number of classical inputs 
in a single run. For instance, starting from the input 
state Y^k^ 1 c k\k)\0), we may obtain the output state 

E <*|*>|/(*)>- ( 4 ) 

fc=0 

Therefore, we have computed the function f(k) for all 
k in a single run. Note that a quantum computer im- 
plements reversible unitary transformations and that we 
need two quantum registers in (4) to compute in a re- 
versible way a generic function /. We emphasize that 
it is not an easy task to extract useful information from 



the output state. The problem is that this information 
is, in a sense, hidden. Any quantum computation ends 
up with a projective measurement in the computational 
basis: we measure the qubit polarization (0 or 1) for all 
the qubits. The output of the measurement process is 
inherently probabilistic and the probabilities of the dif- 
ferent possible outputs are set by the basic postulates 
of quantum mechanics. Given the state (4), we obtain 
outcome \k)\f(k)) with probability \c^\ 2 . Hence, we end 
up with the evaluation of the function f(k) for a single 
k = k, exactly as with a classical computer. Nevertheless, 
there exist quantum algorithms that exploit quantum in- 
terference to efficiently extract useful information. Note 
that in principle one could consider more general, non 
projective measurements. However, this would not mod- 
ify considerations on the efficiency of a given quantum 
algorithm. The reason is that generalized mesurements, 
after the addition of ancillary qubits, can always be rep- 
resented by unitary transformations followed by standard 
projective measurements 1 ' 2 . 

In 1994, Shor discovered a quantum algorithm which 
factorizes large integers exponentially faster than any 
known classical algorithm. It was also shown by Grover 
that the search of an item in an unstructured database 
can be done with a square root speed up over any classical 
algorithm. A third relevant class of quantum algorithms 
is the simulation of physical systems. In particular, the 
power of a quantum computer is intuitive for quantum 
many-body problems. It is indeed well known that the 
simulation of such problems on a classical computer is a 
difficult task as the size of the Hilbert space grows ex- 
ponentially with the number of particles. For instance, 
if we wish to simulate a chain of n spin-i particles, the 
size of the Hilbert space is 2™. Namely, the state of this 
system is determined by 2™ complex numbers. As ob- 
served by Feynman in the 1980's, the growth in memory 
requirement is only linear on a quantum computer, which 
is itself a many-body quantum system. For example, to 
simulate n spin-i particles we only need n qubits. There- 
fore, a quantum computer operating with only a few tens 
of qubits can outperform a classical computer. Of course, 
this is only true if we can find efficient quantum algo- 
rithms to extract useful information from the quantum 
computer. Quite interestingly, it has been shown that a 
quantum computer can be useful not only for the investi- 
gation of the properties of many-body quantum systems, 
but also for the study of the quantum and classical dy- 
namics of complex single-particle systems. As we shall 
see in the next section, basic quantum mechanical prob- 
lems already can be simulated with 6—10 qubits, while 
about 40 qubits would allow one to make computations 
inaccessible to today's supercomputers. We note that 
this figure has to be compared with the more than 1000 
qubits required to the Shor's algorithm to outperform 
classical computations. Moreover, these estimates ne- 
glect the computational overhead needed to correct errors 
and error-correction becomes necessary to obtain reliable 
results by means of a large-scale quantum computer. 
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IV. THE QUANTUM SIMULATOR 

To illustrate the working of a quantum algo- 
rithm by means of a concrete example, let us con- 
sider the quantum-mechanical motion of a particle 
in one dimension (the extension to higher dimen- 
sions is straightforward) 13,14 ' 15 . It is governed by the 
Schrodinger equation 



H^(x,t). 



where the Hamiltonian H is given by 



(5) 



(6) 



The Hamiltonian Ho = — (ft 2 / 2m) d 2 /dx 2 governs the free 
motion of the particle, while V(x) is a one-dimensional 
potential. To solve Eq. (5) on a quantum computer with 
finite resources (a finite number of qubits and a finite se- 
quence of quantum gates), we must first of all discretize 
the continuous variables x and t. If the motion essentially 
takes place inside a finite region, say — d < x < d, we de- 
compose this region into 2™ intervals of length A = 2d/2 n 
and represent these intervals by means of the Hilbcrt 
space of an n-qubit quantum register (this means that the 
discretization step drops exponentially with the number 
of qubits). Hence, the wave function tp(x,t) is approxi- 
mated by 



2"-l 



k=0 



k=0 



where 



x k = -d+ (fc+ f) A, 



(7) 



(8) 



\k) = |fc n -i)|fc n -2) • ■ • \ka) is a state of the computational 
basis of the n-qubit quantum register and 



AT 



E i#»*>*)i 5 



(9) 



fe=0 



is a factor that ensures correct normalization of the wave 
function. It is intuitive that (7) provides a good approx- 
imation to \ip) when the discretization step A is much 
smaller than the shortest length scale relevant for the 
description of the system. 

The Schrodinger equation (5) may be integrated for- 
mally by propagating the initial wave function tp(x, 0) for 
each time-step e as follows: 



1/>(X, t + e) = e -TA H o+V(x)]e ^ f ) _ 



(10) 



If the time-step e is small enough (that is, much smaller 
than the time scales of interest for the dynamics of the 
system), it is possible to write 



-±[H + V(x)]e _ e -± Ho e e -j-V(x)e 



(11) 



Note that this equation, known as the Trotter decom- 
position, is only exact up to terms of order e 2 since the 
operators Ho and V do not commute. The operator on 
the right-hand side of Eq. (11) is still unitary and sim- 
pler than that on the left-hand side. We can now take 
advantage of the fact that the Fourier transform can 
be efficiently performed by a quantum computer. We 
call p the momentum variable conjugate to x, that is, 
—ih(d/dx) — F~ 1 pF, where F is the Fourier transform. 
Therefore, we can write the first operator in the right- 
hand side of (11) as 



F . 



(12) 



In this expression, we pass, by means of the Fourier 
transform F, from the ^-representation to the p- 
representation, in which this operator is diagonal. Then, 
using the inverse Fourier transform F~ l , we return to the 
^-representation, in which the operator exp(— iV(x)e/ft) 
is diagonal. The wave function ip(x,t) at time t = le is 
obtained (up to errors 0(e 2 )) from the initial wave func- 
tion ip(x, 0) by applying I times the unitary operator 



-*(£)« 



Fe -iV(x)e_ 



(13) 



Therefore, simulation of the Schrodinger equation is 
now reduced to the implementation of the Fourier trans- 
form plus diagonal operators of the form 



\x) 



x). 



(14) 



Note that an operator of the form (14) appears 
both in the computation of exp(— iV(x)e/ft) and of 
exp^iHoe/h), when this latter operator is written in 
the p- representation. Any operator of the form (14) can 
be implemented by means of 2™/2 generalized controlled- 
phase shift gates, which apply the transformation Fk to 
the target qubit only when the n — 1 (control) qubits are 
in the state |fc). Here Fk is a single-qubit gate, mapping 
|0) into e i/(2fc)| ) and |1) into e l/ ( 2fc+1 )|l). It is easy to 
check this construction for the three-qubit case shown in 
Fig. 2. F acts only when the first two qubits are in the 
state 1 00) and therefore it sets the phases e 1 -^ ) and e li ^ 
in front of the basis vectors 1 000) and 1 001) , respectively. 
Similarly, F\ acts only when the first two qubits arc in 
the state 1 01) and therefore it sets the phase e*-^ 2 ' and 
e J /(3) m front f the basis vectors |010) and 1 01 1) and so 
on. 

Note that the implementation described in Fig. 2 is 
inefficient, as it scales exponentially with the number of 
qubits. On the other hand, efficient (that its, polynomial 
in n) implementations are possible for most of the cases of 
physical interest. For instance, n 2 two-qubit phase-shift 
gates are sufficient for the harmonic oscillator potential 
V(x) — ^moj 2 x 2 . Using Eqs. (8) and (2), we can write 

the discrctized variable x as a 

constants a = A and 8 = ~ d+A / 2 _ Therefore, x 2 = 

1 rvn. ' 
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-O 



'1 



FIG. 2: A quantum circuit implementing a generic diagonal 
unitary transformation (14), for the case of n = 3 qubits. 
Empty or full circles indicate that the operation on the target 
qubit is active only when the control qubits are set to or 1 , 
respectively. 



a 2 E", il ( k 3 23 + P)^ 1 + 1 3 ) and finall y 



-V(x)e 



IP 

3,1=0 



•* 7 (fc j 2-'+/3)(fc 1 2 , +/3) 



(15) 



where 7 = muj 2 a 2 e / (2h) . This is the product of n 2 phase- 
shift gates, each acting non-trivially (differently from 
identity) only on the qubits j and I. Since the kinetic 
energy Hq is proportional to p 2 , an analogous decompo- 
sition is readily obtained, in the momentum eigenbasis, 
for exp(—iH e/H). Efficient implementations are possi- 
ble for piecewise analytic potentials V(x) but require, in 
general, the use of ancillary qubits. 

Finally, we point out that there is an exponential ad- 
vantage in memory requirements with respect to classical 
computation, While a classical computer needs O(N) bits 
to load the state vector of a system of size N (that is, 
the coefficients ip( x k, t) of its expansion over the compu- 
tational basis), a quantum computer accomplishes the 
same task with just n = log 2 N qubits, namely with 
memory resources only logarithmic in the system size. 
As we have discussed in Sec. II, this is possible thanks to 
entanglement of the qubits in quantum computer. 



A FEW SNAPSHOTS 






FIG. 3: Plots of \tp(x,t)\ 2 with n = 6 qubits (t horizontal 
axis, divided in 40 time steps; x vertical axis, discretized by 
means of a grid of 2 n = 64 points): accelerated particle (top 
left), transmission and reflection through a square barrier (top 
right), harmonic potential, with a squeezed state of initial 
width twice that of a coherent state (bottom left), anharmonic 
potential (bottom right). The initial state ip(x,Q) in all cases 
is a Gaussian wave packet. 



H 



H 



W(t) — W (t) 



FIG. 4: Quantum circuit measuring the squares of the real 
and imaginary parts of the wave function ip(x,t). The line 
with a dash represents a set of n qubits. 



To illustrate the working of the quantum simulator, 
we simulate it on a classical computer, obviously with an 
exponential slowing down with respect to a true quan- 
tum computation. Let us consider a few significant ex- 
amples of single-particle quantum mechanical problems, 
for which plots of \tp(x,t)\ 2 are shown in Fig. 3. In con- 
trast to classical simulations, in quantum computation 
we cannot access the wave function ip(x,t) after a sin- 
gle run up to time t. Indeed, each run is followed by a 
standard projective measurement on the computational 
basis, giving outcome Xk with probability \(k\tp(x,t))\ 2 . 
Therefore, the probability distribution \ip(x, t)\ may be 
reconstructed, up to statistical errors, only if the quan- 
tum simulation is repeated many times. If outcome xt is 
obtained Nj~ times in N runs, we can estimate \ip(xk, t)\ 2 
as Af 2 -^r 7 W1 th Af normalization factor defined in Eq. (9) . 



Fig. 3 exhibits several interesting features of quantum 
mechanics. In the top left plot we can see the spread- 
ing (quadratic in time) of a Gaussian wave packet for 
an accelerated particle. In the top right plot interfer- 
ence fringes appear when a Gaussian packet impinges on 
a square barrier. The bottom left plot shows the width 
oscillations for a squeezed state, back and forth between 
a minimum and a maximum value. Finally, the bottom 
right plot illustrates the motion of a wave packet in a po- 
tential harmonic for x > and anharmonic (V(x) — ax 3 ) 
for x < 0. We can see the deformation and spreading of 
the wave packet when it moves in the anharmonic part 
of the potential. As a result, the initial coherence of the 
state is lost in a couple of oscillation periods. 

As shown in Fig. 4, it is also possible to measure 
the squares of the real and imaginary parts of ip(x, t) 
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by means of a Ramsey type quantum interferometry 
method. For this purpose, we need a single ancillary 
qubit, initially in the state |0). The input state of the 
other n qubits is |0) = |00 • • • 0) and they end up, de- 
pending on the |0) or |1) state of the ancillary qubit, 
in the states \ip(t)) = W(t)\0) or \ip*(t)) = W*(t)\0). 
Note that the operator W(t) = U(t)S includes both the 
state preparation (j-0(O)) = 5|0)) and the time evolution 
(|^(t)} = U(t)\ip(0))). Finally, a Hadamard gate applied 
to the ancillary qubit leads to the (n + l)-qubit output 
state: 



*(*)> = *\o)(\m)+mt))) 



mm) 



(16) 

Hence, the probabilities of obtaining outcomes or 1 
from the measurement of the ancillary qubit and Xk from 
the measurement of the other n qubits are given by 

P (x k ,t) = |(0|(fc|$(i))| 2 = {Rc^xfe,;)]} 2 , 
Pi(x fc ,t) = |<1K*:|*(*)>| 2 = {Im [^(a; fc ,*)]> 2 . 



(17) 
(18) 



From the measurement of {Re [ip(xk, t)]} 2 and 
{Im [i[)(xk, t)]} 2 we can derive the phase velocity from 
the isophase curves of the wave function ip(x,t) — 
\i/)(x,t)\e %<> ( x,t \ corresponding to 6(x,t) = 0,7r (imply- 
ing Im[ip(xk,t)] = 0) and 6(x,t) — ^ir,^ir (implying 
Re [ip(xk,i)] = 0). The phase velocity is just given by the 
slope of the isophase curves. In Fig. 5, {Re [i[)(xk, t)]} 2 
and {Im [i[)(xk, t)]} 2 are shown for the quantum motion 
of a uniformly accelerated particle. It is interesting to 
note that in this example the phase velocity has strong 
local variations: it can change sign and also diverge. 

The huge memory capabilities of the quantum com- 
puter appear in the plots of Fig. 6, where the dynamics 
of the superposition of two counterpropagating Gaussian 
packets in a harmonic potential is considered. When the 
two packets collide interference fringes appear, similarly 
to the double-slit experiment for free particles. These 
fringes are hardly visible with n = 6 qubits but already 
very clear with n = 8 qubits. Since the number of lev- 
els grows exponentially with the number of qubits, the 
discretization step reduces exponentially. Therefore, po- 
sition resolution improves 4 times when moving from 6 
to 8 qubits (of course, there is a lowest resolution limit 
imposed by the Heisenberg uncertainty principle). 

Finally, we point out that the use of the quantum 
Fourier transform in the simulation of the Schrodingcr 
equation automatically imposes periodic boundary con- 
ditions. Such boundary effects become relevant when the 
finite region — d < x < d is not large enough for a cor- 
rect description of the motion. This problem can be seen 
in the two top plots of Fig. 3. In the case of the uni- 
formly accelerated particle, when the wave packet gets 
close to the upper border x = +d, there exists a spurious 
non-negligible probability to find the particle close to the 
lower border x = — d. With regard to the problem of the 
scattering of a particle by a square potential, we observe 
interference fringes between the transmitted and the re- 
flected packets, when both are close to the boundaries. 





FIG. 5: Plots of {Re[ip{x,t)]} 2 and {Im[ip{x,t)]} 2 for the uni- 
formly accelerated particle of Fig. 3 top left. 




FIG. 6: Plots of \tp(x, t) \ for the superposition of two counter- 
propagating wave packets in a harmonic potential, simulated 
with n — 6 (left) and n — 8 (right) qubits. 



In general, it might be interesting for several physical 
problems to consider non-periodic boundary conditions. 
Keeping the quantum Fourier structure of the problem, 
this can be done if suitable fictitious potentials are added. 
For instance, zero boundary conditions may be approx- 
imated if a very high, in practice impenetrable poten- 
tial barrier encloses the region of interest for the motion. 
Twisted boundary conditions [ip(d,t) = e 4< ^?/>(— (i, i)] are 
obtained by adding a suitable Aharonov-Bohm flux <j> to 
the problem. 



VI. FINAL REMARKS 

To asses the feasibility of a quantum computer, 
one should take into account that, in any realis- 
tic implementation, errors due to imperfections in 
the quantum computer hardware or to the unde- 
sired computer-environment coupling unavoidably ap- 
pear. First studies 1 ' 15 have shown that a quantum sim- 
ulator with n = 6 — 10 qubits is quite robust against 
errors. Nowadays, these numbers of qubits are avail- 
able in experiments with liquid state nuclear magnetic 
resonance-based quantum processors and with cold ions 
in a trap 1 . On the other hand, the number of elementary 
quantum gates required for the simulations discussed in 
Figs. 3, 5, and 6 is O(10 4 ), much larger than the se- 
quences of 10 — 100 gates that can be reliably implc- 
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mcnted in present-day laboratory experiments. 

Even though the time when a useful quantum com- 
puter will be realized is uncertain, quantum computation 
is a very promising and fascinating field of investigation 
in physics, mathematics and computer science. Further- 
more, we believe that quantum computation and more 
generally quantum information provide an excellent ap- 



proach to teaching basic quantum mechanics, because 
they deal in an intuitive, appealing and mathematically 
simple way with the main features of the quantum the- 
ory, from the superposition principle to entanglement, 
quantum interference and quantum measurements. 

One of us (G.B.) acknowledges support by MIUR- 
PRIN 2005 (2005025204). 
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